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Abstract. The most general coordinates transformations that allow for the 
exact separation of the kinetic energy operator of a quantum many-body system 
into total center of mass kinetic energy and internal kinetic energy are found 
and discussed. We find i) that the suitable transformations, depending on the 
number of particles, have a certain number of free parameters and this allows for 
the generalization of the Jacobi coordinates to a much larger class of coordinates 
with the same properties and ii) that there is a new, uncommon, additive group 
structure hidden in the transformation matrices that is connected to certain 
geometric properties of the set of coordinates. 
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The kinetic energy operator of a set of particles in quantum mechanics depends 
upon the Laplacian of the coordinates, 

i i 1 

where the summation is extended over all particles i — {1, • • • , n} with position vectors 
ft . The hamiltonian of the system depends on this operator and in the study of the 
internal excitations of a quantum system we are faced with the fundamental problem 
of eliminating the kinetic energy of the center of mass from the total hamiltonian. 
Several transformation of coordinates, techniques and methods exists to accomplish 
this task [TJ [2J, let me just mention the transformation to Jacobi coordinates that 
singles out the center of mass energy [3]. This is a crucial step in atomic, molecular 
and nuclear physics. The general idea behind this is the following conversion 

n n— 1 2=^2 

r = 5>^£i;' + — ( 2) 

i i 

where P — (1/n) Y^i Pn is the linear momentum of the center of mass of the system 
with total mass M = Y^i m i- 

It is generally recognized that there are several sets of Jacobi coordinates that 
accomplish this task [1] and that, not only one can change the labeling of particles, 
but can also change the order in which particles are grouped together, thus generating 
a large number of possible sets of Jacobi coordinates. Then, of course, several people 
have thought to the problem of coordinates transformation for various reasons or 
applications, and they have even included masses (let me cite Ref. 0], where a matrix 
very similar to the ones I derive in the following appears), but to my knowledge 
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the most general conditions under which one can expect center of mass separation 
to occur have not yet been thoroughly worked out. I must therefore apologize to 
all those authors whose researches are not cited in the present paper, but my intent 
is to look at the problem under a wider angle and specific applications, although 
arguably very important, are not referenced. In addition the present analysis reveals 
an uncommon group structure hidden in the transformation matrices that has a 
connection with geometrical properties and might, very likely, open up this subject to 
a re-interpretation. 

In medias res, let's consider coordinate transformations of the type 

x" = Ax (3) 

where the components of the vectors Xj x are connected to each other by 



: ^aijXj . (4) 
3=1 

It is convenient to remember that each of these components might also represent a 
multidimensional vector. 

The relations holding among the corresponding partial derivatives are 

i—l 1 i—l 1 

where we have used 

dx'i d 



dx k dx k ^ 



^ a v x j = X! a ^^i k = a * k ■ ( 6 ) 

3=1 3=1 

We are looking for the conditions to impose on the generic coordinate 
transformation that bring the sum of pure second derivatives of the unprimed 
coordinates into the sum of pure second derivatives of the primed coordinates without 
introducing mixed derivatives. In other words for each k = 1, ■ • -n, if 

Ql / n Q \ / » Q \ " Q2 

^ = \JL a *Mj (J>*azJ = ? a ^ ( 7 ) 

» x i—l L/ j — 1 3 — 1 1 J 

is the transformation of the second derivatives, the n(n — l)/2 conditions to forbid the 
presence of mixed derivatives are 

n 

^2a lk a jk = V{i,j} /i ^ j = {1, • • • , n} . (8) 
fc=i 

These n(n — 1)/2 equations form a system with n 2 unknowns, hence n(n + 1)/2 matrix 
elements might arbitrarily be chosen. In order to insure that the last of the coordinates 
always coincides with the position of the center of mass, the elements in the last row 
must take the form a„i = 1 /n, V« = 1, • • • , n, and hence upon simplification the system 



is replaced by 



ELi a ^ a ife=° 1% ^ j = {l,---,n- 1} 

(9) 

J]fe=i a ifc = i = {1, • • ■ ,n - 1} 



where the first line represents a subsystem of {n — 2)(n — l)/2 conditions of the form 
([8]) with indexes i,j restricted to avoid n and the second line represents n — 1 simpler 
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conditions. In general the system §§§ has n 2 — n unknowns, of which (n — — 2)/2 
and (n — 1) are set and the remaining matrix elements can be chosen arbitrarily. 
Therefore there are entire classes of transformations that comply with the requirement 
(J5J. Oftentimes it is useful to have that at least some of the vectors to correspond 
with relative position vectors and this implies transformation matrices with (in some 
rows) only two elements different from and equal to 1 and —1 respectively. 
Example: for two particles, n = 2, we have the matrix 

( an ° 12 ^ (10) 
The condition ((SJ) takes (twice) the form 

011021 + 012022 = 

that, remembering that we wish to keep 021 = 022 = 1/2, has only one solution 
an = -an . 

Leaving aside the trivial an = at2 = 0, that is unacceptable to our scopes, 
we have the freedom to choose the value of one of the matrix elements. Note that 
with an = 1 it is exactly the coordinate transformation that takes the position of 
two particles into the relative coordinate plus the coordinate of the center of mass, 
depicted in Fig. ([1]) . Any other choice would have simply altered the modulus or sign 
of the internal coordinate x' 2 without changing its direction. 




Figure 1. Transformation discussed in the text for 2 particles (black dots). 
Hollow circles represent the external reference point and double circle represents 
the particles' center of mass. Dotted lines are rails (see text). 



The most general matrix we can use for 2 particles is 

*<°>=(l/2 1/2) (u) 

with freedom on a. The transformation of coordinates of Fig. {I) corresponds to 
,4.2(1). Note that these matrices form an additive group under the law of composition 
obtained by the operation of standard matrix addition and subsequent division by 
two. Although there is nothing exceptionally new about this composition law (it is 
essentially matrix addition), I will indicate it with the -K- symbol, for the sake of 
brevity. Instead these matrices are not a group under usual matrix multiplication. 
In other words, we don't try to apply ^2(0) to ^(a) and then to the column 
vector (xi,X2)- by doing so one ends up with something completely useless. But 
rather apply the composition A^ia) +-^^2(6) to the column vector: this will give 
again a transformation of the same type. The meaning of this group is that all the 
infinite possible transformations that preserve the pure second derivatives and separate 
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the second derivative with respect to the center of mass coordinate are essentially 
equivalent, and therefore the parameter a can be set to 1 without loss of information. 
The entire class can be seen geometrically as follows: in Fig. |T]) draw the two lines 
connecting the external reference point with 1 and 2 and call them rails; shift the 
vector x[ along the rails, keeping it parallel to the original and adjusting the length 
in order to maintain its extremes on their rails. The infinite set of vectors that are 
obtained by this parallel expansion/contraction procedure maps one to one to the real 
numbers a. 

Example: for three particles, n = 3, the matrix is 
an ai2 ai3 \ 

021 022 023 (12) 
^31 «32 033 / 

and the conditions are 

OuOai + a 12 a 22 + a 13 a 23 = 

ono 3 i + a i2 a 32 + ai 3 a 33 = (13) 

021031 + 022032 + 023033 = 

that with CI31 = a32 = 033 = 1/3 can be simplified to a system of three equations with 
six unknowns 

a n a 2 i + a 12 a 22 + a 13 a 23 = 

an+ai2+ai3 =0 (14) 

021 + 022 + 023 =0 

where two of the three equations have the form suggested in (J9j) . Therefore one has 
the freedom to choose three elements and deduce the others according to (|14p . 




Figure 2. Transformation discussed in the text for 3 particles (black dots). 
Hollow circles represent the external reference point and double circles represent 
the particles' total and partial centers of mass. 



Using, for instance, an = 1, a\ 2 = —1 and 022 — 1/2 one gets for x! a possible set 
of Jacobi coordinates containing the position of the total center of mass, the relative 
coordinate between particle 1 and 2 and the relative coordinate between particle 3 
and the minor center of mass of particles 1 and 2, as shown in Fig. @. 

The set of equations (fT4"| is important, because it tells at a glance that, for 
example, certain coordinates systems don't have the right property of the separation 
of the center of mass. Consider the coordinates as in Fig. ([3]), sometimes called 
necklace or sequential coordinates Q. Using an = 1, a\ 2 = —1 one gets 013 = 0, 

I A folklore note: it is precisely the fact that i was considering the necklace coordinates in connection 
with many-body nuclear physics models that brought my attention to to the fact that this problem 
was not thoroughly known and to the formulation of the general conditions under which the separation 
is possible. 
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0-21 = &22 and therefore 2a 2 2 + ^23 = 0- While the choice of the preceding example 
works well, the choice et 2 2 = 1 and a 2 3 = — 1 does not satisfy this last equation and 
therefore does not permit the separation of the center of mass as described by @. 




Figure 3. A coordinate system that does not obey the set of equations l|14|l and 
thus does not allow for the separation of the center of mass kinetic energy. 



The most general matrix for 3 particles is 

(a b —a — b 

c i " : " 

1/3 1/3 1/3 

with arbitrary a, b, c. These matrices are not in general a group as before, but 
remarkably they can still form an additive group under the operation -K- by restricting 
some of the parameters. This can be accomplished only in two ways: either b = ka 
or a — k! c — 26, where k and k! are constants. By imposing these two conditions the 
resulting transformation matrices are respectively 

i(l + k) \ 

(16) 





( a 
V 1/3 


ka 

r 2+fe 
°l+2fc 

1/3 


-a(l + k) 

L l+2fc 

1/3 




/ k'c- 


2b b 


b- 




A«(6,c)* = c -2c+f c-f I (17) 

V 1/3 1/3 

where the superscripts in parenthesis just denote an arbitrary labeling. The difference 
between these two matrix groups lies in the fact that one has two arbitrary parameters 
in the same column, and the other in different columns and, upon reshuffling of 
the columns, these are the only two possible ways. It can be noticed that the 
transformation to Jacobi coordinates of Fig. @ is an element of both groups that 
corresponds to 

4,(1, -1. 1/2) - 4 1} (1, l/2)-i = 4V1, l/2)-2 ■ (18) 
An alternative set of Jacobi coordinates is given by ^(0, 1, —1). Notably this is not 
an element of A3 , but it is an element of A 3 2 '(l, — 1)_2- However, by noticing that its 
matrix corresponds to the previous Jacobi coordinates matrix upon reordering of the 
columns, one can see that this fact is simply explained by the arbitrariness in choosing 
the free matrix elements in (|15p : we have filled the upper left triangle with arbitrary 
parameters and determined the rest, but one could as well have made other choices. 
Summarizing, we must keep in mind that there are two levels of arbitrariness in the 
procedure of calculation of the matrices: i) it is arbitrary to choose which elements 
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are parameters and which are functions of these parameters and ii) the values of the 
parameters themselves are arbitrary. 

Several choices of parameters a, b, c can be traced back to a "rail projection" of 
a set of Jacobi coordinates, as, for instance, ^(l, 1, 1): the first coordinate coincides 
with x[ of Fig. ([2]), the third is the center of mass coordinate x' 3 as in all other cases, 
but the second is just an expansion of x' 2 of Fig. ([2|) along the two rails connecting 
the reference point with 3 and with the center of mass of particles 1 and 2. This 
is analogous to the case already discussed for 2 particles. Apart from these cases, 
there are others that, despite their strange look, are equally good at separating out 
exactly the kinetic energy of the center of mass. Take for example ^(1/2, 1,-1) 
that gives the coordinate transformation depicted in Fig. (j4]). Albeit this might be 




Figure 4. A rather weird coordinate system corresponding to As(l/2, 1, —1) that 
nevertheless allows for the separation of the center of mass kinetic energy. 



useless to the solution of the many-body problem, it is a viable alternative to Jacobi 
coordinates. This example shows clearly that, when n > 2, the group structure has a 
geometrical interpretation that is wider than the parallel expansion/contraction along 
the rails described above, but several subsets of all possible transformations still have 
this interpretation. In particular for the matrix (|16p it is possible to find fixed rails, 
because a and c amount to a multiplicative constant for x[ and x' 2 respectively, but in 
this case, by changing k, one is not only shifting the vector's tips along the rails, but 
it's also changing its direction accordingly. The interpretation of (fT"7|) is instead less 
obvious, because two parameters enter at the same time into several matrix elements. 

For n particles we return to the system The corresponding transformation 
matrix that maps the initial set of coordinates into a set of coordinates that allows 
the exact separation of the center of mass kinetic energy can be schematically divided 
in the following way: 



(19) 



\ X/n 



cii; 



where the n matrix elements of the last row are all equal to 1/n. Conditions (j9]-2) is 
used to determine, for example, the remaining n — 1 elements of the last column (one 
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could have equivalently chosen another column): 

n-l 

Pin = - ^ <Hk i = {1, ■■■,n- 1} . (20) 

k=l 

Conditions (J5J-1) are used to set the (n — l)(n — 2)/2 matrix elements of the lower 
triangle as functions of the arbitrary p = n(n — l)/2 matrix elements of the upper 
triangle, starting from the pair = {1,2} and proceeding row by row. Matrices 

P9[) most certainly don't form a group under 44-, but it is still possible that, as in the 
case n = 3 they would do so by restricting the number of parameters under certain 
special conditions. 

Let us give an example of the usefulness of the conditions ©. With reference to 
Fig. ([5]), suppose that we wish to use the internal coordinates of the two subsystems 
formed by particle 1 — 2 and 3 — 4 and in addition we want also the right property of 
separating the total kinetic energy center of mass (into the new coordinate x^). 




Figure 5. Coordinate system corresponding to A±(l, — 1, 0, 1/2, 1/2, 0). The rails 
are indicated as dotted lines. This system can be used to model scattering of two 
deuterons for example. 

By filling up the transformation matrix for four particles with the proper values 
for the first, third and fourth row, we get 



/ 1 


-1 





° ^ 


a 


a 


—a 


—a 








1 


-1 


V 1/4 


1/4 


1/4 


1/4 J 



(21) 



therefore the entire class of possible transformations in this particular case reduces 
to all the vectors with direction parallel to that the center of mass of 3 — 4 to the 
center of mass of 1 — 2, or, in other words all expansions/contractions along the rails 
shown in figure as dotted lines. With the choice a = 1/2 we have the simplest of these 
vectors, i.e. x' 2 of Fig. ([5]). This coordinates are suitable to describe, for instance, the 
scattering of two deuterons, where the total wave function can be expanded in terms of 
products like ^12 (5a) 034 (#3) ^(^2)) wne re <fi are internal wavefunctions that describe 
separately the subsystems 1 — 2 and 3 — 4 and ip is the relative motion wavefunction 
between the two subsystem. Several other systems might be treated with coordinates 
of this type: cluster states in light nuclei, tetratomic molecules and so on. 

In conclusion, we have discussed the most general conditions to impose on 
coordinates transformations that allow for the exact separation of the total center 
of mass kinetic energy, obtaining the set of equations ©. These equations show that 
there are classes of coordinates systems that can be used in place of usual Jacobi 
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coordinates, achieving the same goal of elimination of the center of mass energy. The 
application to small numbers of particles has been discussed and we have found that 
a curious additive group structure is hidden in these transformation matrices. We 
have analyzed it finding a connection with a simple geometrical interpretation. We 
believe that these results can be of some importance to quantum many-body models 
based on the solution of the non-relativistic Schrodinger equation and that the new 
group structure might reveal further hidden symmetries and analogies that might help 
advancing this field. 
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